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• The He2-SF6 trimers, in their different He isotopic combinations, are studied both in the frame- 

' work of the correlated Jastrow approach and of the Correlated Hyperspherical Harmonics expansion 

method. The energetics and structure of the He-SFe dimers are analyzed, and the existence of a 
characteristic rotational band in the excitation spectrum is discussed, as well as the isotopic differ- 
ences. The binding energies and the spatial properties of the trimers, in their ground and lowest 
lying excited states, obtained by the Jastrow ansatz are in excellent agreement with the results of 
O . the converged CHH expansion. The introduction of the He-He correlation makes all trimers bound 

Jh ' by largely suppressing the short range He-He repulsion. The structural properties of the trimers are 

S \ qualitatively explained in terms of the shape of the interactions, Pauli principle and masses of the 

. constituents. 
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I. INTRODUCTION 



Helium systems are dominated by quantum effects and remain liquid down to zero temperature. This is a conse- 
quence of both the small atomic mass and the weak atom-atom interaction, which is the weakest among the rare gas 
atoms. Helium clusters remain liquid under all conditions of formation and are very weakly bound systems. 
, Small helium clusters have been detected by diffraction of a helium nozzle beam by a transmission grating, ?^ Using a 
■ grating of 200 nm period, conclusive evidence of the existence of the dimer ''He2 has been established.^^ The existence 
T-H of ^He2 was previously reported^ using electron impact ionization techniques. Diffraction experiments from a 100 
0^ nm period grating has lead to the determination of a molecular bond length of 54 ± 4 A, out of which a binding 
energy of 1.1 0.3 - 0.2 mK has been deduced.-* This energy is in good agreement with the values obtained by direct 
• integration of the Schrodinger equation using modern He-He interactionsi^ Theoretical calculations indicate that any 
number of "^He atoms form a self-bound system. In contrast, a substantially larger number of ^He atoms is necessary 



for self-binding, as a consequence of the smaller "^He atomic mass and its fermionic nature. The required minimum 
. ^ number has been estimated to be 29 atoms, using a density functional plus configuration interaction techniques to 
£^ , solve the many-body problem^ or 34-35 atoms, using an accurate variational wave functionL^ with the HFD-B(HE) 
^H ^' Aziz interactioufS It is worth recalling that a theoretical description of pure ^He clusters, either based on ab initio 

calculations or employing Green Function or Diffusion Monte Carlo techniques, is still missing. 
• • . Doping a helium cluster with atomic or molecular impurities constitutes a useful probe of the structural and 
. !^ ' energetic properties of the cluster itself. It has been proved that rare gases and closed-shell molecules as HF, OCS 
^ or SFg are located in the bulk of the cluster liS Doped '^He clusters have been extensively studied in the past by a 
^ variety of methods, ranging from Diffusion and Path Integral Quantum Monte Carlo methods to two-fluid models. 
. 5^ , A comprehensive view of this subject is found in Ref. illi giving account on a microscopic basis of the free rotation 
of a heavy molecule in a "^He nanodroplet, consistently with the occurring bosonic superfluidity at the attained 
temperatures. This phenomenon is not expected to take place in doped '^He clusters, because of the fermionic nature 
of the atoms, unless very low temperatures (below a few mK) are reached. Even more interesting, from both the 
theoretical and experimental points of view, are mixed doped clusters of '^He and ^He atoms. The status of the theory 
in these last cases is far behind that in the *He droplets. The most updated studies of ^He and mixed *He-'^He 
clusters, either pure or doped, employ a finite-range density functional theory^Sii^ 

In this work we study the properties of a trimer formed by two helium atoms plus a heavy dopant. The dopant 
molecule behaves as an attractive center binding a certain number of otherwise unbound "^He atoms. This fact has 
been used in Ref. to set an analogy between electrons bound by an atomic nucleus and ^He atoms bound by a 
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dopant species. Systems formed by two ^He atoms plus a molecule have been studied employing the usual quantum 
chemistry machinery. 

Helium droplets doped with the SFg molecule have been widely investigated, also in view of the fact that the 
interaction He-SFg is well establishedii^ti^ In this work we use the spherically averaged interaction of Taylor and 
Hurly^^ between the helium atom and the SFg molecule, and the Aziz HFD-B(HE) helium-helium interactioniSi We 
first perform a variational study of the '^He2-SF6, ■*He2-SF6 and ^He-'^He-SFg trimers using a Jastrow correlated wave 
function. Then the Correlated Hyperspherical Harmonics (CHH) expansion methodii is employed and its outcomes 
are used as benchmarks for the variational calculations. The CHH expansion has shown to be a powerful technique 
to study three-and four-body strongly interacting systems. In light atomic nuclei its accuracy is comparable with 
(and in some cases even better than) other popular approaches, as the Faddeev, Faddeev-Yakubowsky and Quantum 
Monte Carlo ones»i^ Besides accurately studying the ground and first excited states of the trimers, comparing the 
variational and CHH results may provide essential clues for the construction of a reliable variational wave function to 
be used in heavier doped nanodroplets. 

The plan of this paper is the following. In Sect. II we study the dimcrs formed by a single helium atom and the SFe 
molecule and enlighten some aspects of their excitation spectrum. In Sect. Ill we consider the trimers by the Jastrow 
variational and CHH approaches. Results for the energetics and structure of the trimers are given and discussed in 
Section IV. Finally, Sect. V provides the conclusions and the future perspectives of this work. 

II. THE HE-SFe DIMER 

Prior to the study of trimers made by two helium atoms and a dopant species, it is convenient to analize the dimer 
in some details. To fix the notation for later discussions, we write here the Schrodinger equation for the relative 
motion of a helium atom and a dopant, D, 

hQ(t>ni{r) = e„f(/)„^(r) (1) 

with 

ho ^-7; V2 + 1/He-D(r) (2) 

Znia 

where nia (with a = 3,4) is the reduced mass of the "He-D pair, and VHo-D(r) is the helium-dopant interaction, 
being r the relative coordinate. We have numerically solved this equation for the dopant SFg, using the spherically 
averaged interaction determined in Ref. A set of energies and orbitals characterized by the quantum numbers {ni) 
is thus obtained. 

The He-SFg interaction has an attractive well strong enough to sustain 12 and 15 bound states for isotopes '^He 
and ''He, respectively. In Table some observables of the dimers are displayed. The calculations have been done in 
the limit of infinite mass of the SFg molecule. The most striking feature of the energy spectra is that the first 7 (9) 
levels correspond to nodeless n = states. Notice that for each isotope the expectation values (R) and y/{B?) are 
not very different, neither for a given ^-state nor for different ^ values. The expectation value {V) neither varies too 
much, increasing by ~ 10% in going from the ^=0 to the ^=8,9 states. These results are an indication that the wave 
functions show a well-defined peak in nearly the same region, as a consequence of the characteristics of the He-SFg 
interaction. 

To first order in the mass ratio, the correction to the infinite dopant mass approximation modifies the kinetic energy 
as T*^"^ — T'j|f2,oo(-'^ + fna/M), with M the dopant mass. Accordingly, the finite mass system results less bound by 
^{a) _ ^ j[,fy fpjjg niass ratios iria/M are 0.0207 and 0.0274 for -^He and ''He, respectively, considering the 

■^^S isotope. So, the binding energy corrections are A^^ = 0.245 K and A^'' — 0.300 K, both being about 1% of the 
total energy. An exact finite mass calculation confirms this estimate, providing A^-* = 0.242 K and A^^ = 0.296 K. 

The He-SFg interaction is displayed in Fig. It is strongly repulsive at distances shorter than ~ 3.8 A, and 
attractive beyond. The attractive part is mostly concentrated in a narrow region around 4.2 A, immediately after 
the repulsive core. This implies that the He atom locates relatively far away from the potential origin, so that the 
centrifugal term entering the Schrodinger equation can be considered as a perturbation. Two consequences can be 
deduced from this observation. First, the radial distributions are very peaked in the same narrow region, independently 
on the value of the angular momentum £, as it is shown in Fig.^for three states, corresponding to ^ = and 4 for both 
isotopes, and the nodeless bound states with the highest excitation energy, namely ^ = 8 for '^He and 9 for ^He. The 
distributions are concentrated in the same region, and those corresponding to £ = and 4 are barely distinguishible. 
Second, the excitation energies are closely proportional to £{£ + 1), i.e. they follow a rotational pattern. Figure [21 
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depicts the differences {eoi — cqq) in functions of £{£ + 1) (squares for ^He-SFg and stars for ^He-SFg), together with 
the hnear fits to the energy differences. The slopes provide the n = rotational constants, Cq'^'"'-'. From the the fits 
we obtain Cq^'' = 0.376 K and Cg'*'' = 0.294 K, in good agreement with the rough estimate /2m{R^), where (i?^) is 
the mean square Dopant-Helium distance. A similar behavior is found for the n — \ excited states, whose rotational 
constants are cf^ = 0.194 K and c'^'^ = 0.179 K. The ratio of the rotational constants of the dimers with either 
isotope are of course in the inverse ratio of masses. Notice that the larger mass of the ^Hc atom translates into a 
larger binding. Both the decrease of the kinetic energy and the increase of the attraction contribute to the increment 
of the binding energy for the ''He-SFg dimcr. The stronger localization of the ^He atom is visualized by the peak of 
the radial probability density shown in Fig. ^ which is slightly higher for the ''He atom. 

It is worth stressing that the observed rotational spectrum is a direct consequence of the shape of the He-SFe 
interaction, whose attractive well, having a depth of « —57 K, allows for twelve or fifteen bound states. Taking 
the lighter Ne as a dopant, we have found that only three bound states exist, being the depth of the attractive well 
« —20 K. Moreover, the binding energies are very small, and the probability distributions are extended over a large 
region. 

III. THE HEz-SFfi TRIMERS: THEORY 

In the limit of infinite mass of the dopant molecule, hence considered as a fixed center, the Hamiltonian of the 
He2-SF6 trimers is 

2 *2 2 

H{1, 2) = - ^ + ^ Vb-Hc(r.) + T4ic-Hc(ri2) (3) 

i—l i—1 

with tti = 3, 4. VD-Ho(r) and yHc-Hc(r) are the dopant-helium and helium-helium interaction potentials, respectively. 
Vi is the coordinate of the i-th helium atom with respect to the central molecule and Vi2 is the helium-helium relative 
coordinate. 

The ground and excited states properties of the trimers can be obtained either by an exact solution of the Schrodinger 
equation, 

i/(l,2)vl/^(l,2) = i?^vl,^(l,2) (4) 

or by some approximate estimate of their wave functions. Here 7 labels the generic trimer state, whose wave function 
is vl/^(l,2). 

The Schrodinger equation for clusters of *He atoms may be exactly solved for the ground state by quantum Monte 
Carlo (QMC) methods ji^tSfiiSi Other approaches, as Variational Monte CarloSSiSi24 (VMC) with Jastrow correlated 
wave functions or density functional theorie3SS*2& (DFT), provide a less accurate description of the clusters. However, 
they are generally more flexible than QMC. 

The presence of more than two '^He atoms makes the exact solution of the Schrodinger equation much more 
difficult, because of the notorious sign problemSl associated to their fermionic nature. As a consequence, only DFT 
based studies of doped '^He clusters are available in literatureji^*^ The most updated study of the ^He2-SF6 trimer 
has been done within the Hartree-Fock (HF) approximation, not considering the strong He-He correlations induced 
by ^Hc-Hc'^ 

In this Section we will first present a variational approach based on a Jastrow (J) correlated wave function. Then we 
will apply the Correlated Harmonical Hyperspherical (CHH) expansion methodic to further improve the description 
of the doped trimers. 

A. Variational approach 

The Jastrow correlated wave function of the trimer for the 7-state is given by: 

vl/:j(l,2) = $^(l,2)/j(ri2) (5) 

where $7(1, 2) is an independent particle (IP) wave function of the two helium atoms in the dopant field having the 
same set, 7, of quantum numbers. The correlation function between the two atoms, fj{r), is assumed to depend only 
on the interatomic distance and takes into account the modification to the IP wave function mainly due to the He-He 
interaction. The optimal fj{r) is variationally fixed by minimizing the total energy of the state. $-,.(1,2) is built as 
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an appropriate combination of the dimer He-SFg wave functions. For instance, using the (LS)'" notation (L being the 
IP orbital angular momentum, S the total spin and vr the parity of the helium pair), the (00)''" IP wave function for 
the ^He2-SF6 trimer is taken as: 

$(oo)+(l,2) - So(l,2)0i,(ri)(/)i,(r2) (6) 

where So(l,2) is the spin-singlet wave function of the '^He-'^He pair and 0is(r) is the Is (n=0 and £=0) solution of 
the •^He-SFe dimer Schrodinger equation Q. 

In an analogous way, the (10)~ and (11)^ IP wave functions are: 



and 



$(10)- (1,2) = So(l,2)-i=[0i,(ri)(/.ip(r2) + (/.ip(ri)(/.i,(r2)] 



$(11)- (1,2) = Si(l,2)-^[0i,(ri)(/.ip(r2)-(/.ip(ri)(/.i,(r2)] 



(7) 



(8) 



where Si (1,2) is the spin-triplet pair wave function. 
The total Hamiltonian ||2J) can be written as: 

Hil,2) ^ hoil) + ho{2) + VHe-He{ri2) ^ Ho{l,2) +VHe-He{ri2) (9) 

where ho is the Hamiltonian Q of the dimer, and 

i/o(l,2)$^(l,2) = ( )<i>^(l,2) (10) 

where the two sets of dimer quantum numbers, 71 and 72, are those taken to build up the total 7-trimer state. 
The total energy of the trimer in the 7-state is 

Since [Vd-Hg, /,/] = 0, the energy results to be: 

p _ _^ ,^ {^,\y^fJirl2)■y^fJ{rl2)\<^,} , ($7 1 // (^^12 )^ge-ge (ri2 )/j (ri2 ) | $7) ^ 

This equation will be used to estimate the variational energy of the trimer and to optimize the choice of the correlation 
factor. 



B. CHH approach 

In order to implement the CHH method for a system of three atoms of masses rrii, in positions r,;, it is convenient 
to introduce the three sets (xi,yi) of Jacobi coordinates: 



_ I mi(mj+mk) / _ mjrj+mkrk \ 



In the fixed-center limit (ma = oo) the position ra coincides with the molecular center of mass and, for two equal 
mass atoms (mi — 1712 = m), the Jacobi coordinates, after dividing by m^/^, can be reduced to: 



xi = r2 X2 = -ri 



= . (14) 

y'='' 1^2 = "^ y3 = -^/^(rl + r2) 



The total wave function, ^(^5)7^(1, 2), can be expressed as a sum of three Faddeev-like amplitudesji^ each of wich 
explicitely depends upon a different Jacobi set: 

^^(^5)^(1, 2) = V(%'(^i'yi) + V';LS)'(^2,y2) + V'(%,(x3,y3) . (15) 
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The amplitudes are then expanded into channels, labelled by the partial angular momenta, i^.i and associated 
with x,j and y^, respectively: 

^jl^5)„(x„yO = Ss(l,2) ^ <^Z,eJx.,y^)[Yi^^^{x,)Y,^Jy.)]LM, (16) 

where Yg are ordinary spherical harmonics, and <i>^*'' . ^ . is a two-dimensional function depending upon the moduli of 
the Jacobi vectors. As a result, the parity of the state is given by £x.i + even (odd) for positive (negative) parity. 

Following Ref. [2^ each amplitude , ^ is expanded in terms of a correlated hyperspherical harmonics basis set. 
After introducing the hyperspherical coordinates, {p,(f)i), associated with the Jacobi set: 

2 2|22|22|2 /'1'7^ 

Xi = pcos(/)i, (18) 
yi = psint^i, (19) 

the CHH basis elements, having quantum numbers {LSy and corresponding to the set of Jacobi coordinates labelled 
by i, are defined as follows: 

|m,fc,4,^,^„,^;i) = Ss(l,2)Fj(ri,r2,ri2)p^-'+^-L(^)(z)exp(-^) 

X ^^^Pt"''''{<l^^)[Yi^Ai^)Y,,^yi)]LM, (20) 

where ^^•'P^'"'^" (</>) — Ni_^^i^^k{cos(j)Y^ {sm(j)Yy P^^^^^'^'^"^^^'^ {cos{24>)) is a normalized hypherspherical polynomial, 
p^„+i/2/x-i-i/2 ^ Jacobi polynomial, are associated Laguerre polynomials, and z = /5p, with j3 a non-linear 
variational parameter. 

The other ingredient in the CHH basis elements is the correlation factor, Fj (ri , r2 , ri2 ) . Due to the strongly repulsive 
core of the interatomic potentials, it is convenient to take Fj as a product of pairwise Jastrow-like correlation functions, 

Fj{r^,r2,r^2) ^ g{ri)g{r2)fj{ri2) . (21) 

The correlation functions mainly describe the short-range behavior of the wave function as the two helium atoms 
are close to each other or to the dopant. The polynomial part of the CHH expansion is expected to reproduce the 
mid- and long-range configurations. Therefore, it is preferable to choose correlation functions which behave in the 
outer portion of the Hilbert space as smoothly as possible, in order to avoid potentially deleterious biasing of the 
polynomial expansion. 

Suitable correlation functions to be used in (|21|) are obtained by solving the two-body Schrodinger equation: 

4-2 \ 

_ V2 + l/*(r) + Ae-«'- /i(r)=0, (22) 

where /i is the reduced mass of the considered two-body system, and the pseudo-potential, Ae^^'', is introduced to 
adjust the asymptotic behavior of the correlation function in such a way that h{r —f 00) — > 1. The parameters ^ 
and A are optimized for each of the different cases, SFg-He, "'He-"'He and ^He-'^He. For the helium-dopant correlation 
{h — g), we take /x = mnc, and V* is a modified SFg-He potential. In fact, in order to build a nodeless correlation 
function, it is necessary to reduce the attractive part of the helium-dopant interaction. The repulsive part, on the 
other hand, is kept unaltered. In the '*He-''He case (/i = /j, /i = the correlation function has been obtained 

as in Ref. 123. The '^He-'^He pair does not support a bound state, so the correlation function has been obtained simply 
as the solution of the zero-energy Schrodinger equation (V^* = Vne-He and A = in Ea. l^ . 

In order to work with basis elements of defined symmetry under the permutation operator (according to the 
Pauli principe) we take proper symmetric or antisymmetric combinations of the basis elements with i = 1 and 2: 

|m, fc, 4,1, v) = k, 4a, ly,i] 1) + (-1)"!™, fc, 4,1, lys] 2) , (23) 

where v (= s, a) labels symmetric or antisymmetric states, respectively, and <7=0, 1 is choosen according to the values 
of 4 and S. 

The expansion of the total wave function in terms of CHH states with well defined permutation symmetry results 



(24) 
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The sum over the partial angular momenta ix,iy^ although constrained by the values of the total angular momentum 
L, the parity and by symmetry considerations, runs over an infinite number of channels. However, in practice only the 
lower channels are included, since the higher the angular momentum the less the contribution to the wave function. 
Moreover, due to the presence of both the correlation factors Fj and the amplitude expansion, an infinite number of 
channels is automatically included, though in a non-flexible way. 

In the description of the [LSY — (00)+, (10)~, (11)^ states for ^He2-SF6, we have retained only the lowest angular 



momentum channels, that is. 



for (00)+, 4 



and I 



ix,3 = 1 and (.y^i = iy^3 = for (00) . The CHH wave functions for these three states are: 
vI/(oo)+(l,2) = So(l,2)^^je-^/2 J ^ Afc„,„„L(^)(z)P,f ■i/^(cos20i) 

I ka ,m 

ktj.nii, 

vI/(io)-(l,2) -So(l,2)^^,,e--/2 J J2 Ak^.rnAnli^) siu 0iPff ' (cos 2^i ) Y^^^ (yi ) 



1 for (10)-, 



(25) 



+pcos,^i Pff ''/'(- cos 2</)i)y*^(£i)J + J2 i?fc.,™.4'2(^)P«in,^3Pfcr''^'(cos2</)3)y*^(y3 



*(n)-(l,2)=Si(l,2)P,,e-^/2<^ ^ Ak^^^^Li^liz) 



F^'/'^^/'(cos20i)pcos0iyi^^(£i) 



ka jf^t'a 



-psin0i P^Y'' (-cos20i)r/^(yi)J + ^ i?fc„„,L(^)(z)pcos</>3P,Y' (cos203)r*^(i3) 

fci,,m(, 



(26) 



(27) 



where fc^ > and fcf, > 1. However, only even values of ka are allowed in 5'(oo)+- In ^He2-SF6 we just study ^"(00)+ 
and \l/(io)- states, since only S = combinations are possible. 

In Eqs. I|25l - I27|l the linear coefficients, {^fe^.m^} and {Bkb,mi}, are unknown quantities to be determined. The 
implementation of the variational principle for linear variational parameters leads to a generalized eigenvalues problem 
whose solutions, E^^ are upper bounds to the true energy eigenvalues of the three-body Schrodinger equation. It is 
possible to improve the estimates of the energies by including a larger number of polynomials and channels in the 
CHH basis. If N is the dimension of the basis set, the estimates will monotonically converge from above to the 
exact eigenvalues as N is increased. The pattern of convergence for the (00)+ two lowest lying state of ^He2-SF6 is 
shown in Table ^] An optimum choice of the non-linear parameter /3 has been adopted to improve the convergence 
rate. 



IV. RESULTS 



Table IIIII collects the energies for the He2-SF6 trimers in the uncorrelated, Jastrow and CHH approaches. The 
correlation function is set equal to unity in the uncorrelated calculations, whereas it has been choosen of the McMillanSS 
form, 



fj{r) = exp 



1 Act 

2 [V 



(28) 



in the variational case. Here, (7=2.556 A and b is the only non-linear variational parameter. This type of correlation has 
been widely adopted in variational studies of liquid helium since it provides an excellent description of the short-range 
properties of the correlated wave function. In fact, it gives the exact short-range behavior for a 12-6 Lennard- Jones 
atom-atom interaction. The correlation operator adopted in the CHH expansion has been described in the previous 
section. In the CHH case we use different correlation functions since a non-linear parameter, is already present in 
the basis functions (see Eg. I20|l . Therefore, employing the McMillan form would imply a two non-linear parameters 
minimization. However, it has been checked that the use of the McMillan correlation in the CHH expansion produces 
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binding energies within 0.01% of those given in Table IIIII showing that the converged results are to large extent 
independent of the correlation function, provided the short range behavior is adequately described. In Table ITm we 
also show the kinetic (T) and potential {V) contributions to the energy, separating the latter in its dopant-helium 
(D-He) and helium- helium (He-He) parts. 

The uncorrelated trimers are unbound in all states, with the exception of the spatially antisymmetric (11)^ for 
'^He2-SF6. The orbital antisymmetry reduces the probability of configurations having the two '^He atoms close to each 
other. Hence, the contribution of the strong He-He repulsion at short distances is drastically suppressed. In Ref. ^3 
an expansion of the single-particle helium s- and p-orbitals in a finite set of gaussian basis functions centered at the 
dopant provided an energy of -31.36 K for the (11)^ state in ^He2-SF6. This energy is higher than our uncorrelated 
estimate, pointing to a lack of convergence in the Hartree-Fock result in that reference. 

The introduction of the Jastrow correlation bounds all trimers, as it suppresses the short-range helium-helium 
repulsion. The L=0, positive parity states are the lowest lying ones, the other states having small excitation energies, 
lower than 1 K. The values of the variational parameter giving the minimum energies are also reported in the Table. 
The 6- values for the spatially symmetric trimers are close to those found in the Jastrow correlated studies of bosonic 
liquid ^He. As in fermionic liquid '^He, b is smaller for the spatially antisymmetric (11)~ state, since both the 
correlation and the Pauli principle concur in depleting the Vnc-Hcif) repulsion. 

The converged CHH expansion provides slightly more binding (at most about -0.1 K) to the trimers. This fact 
is a strong indication of the high efficiency of the simple Jastrow correlated wave function in these systems. The 
kinetic and dopant-helium potential energies do not vary much in going from the uncorrelated to the variational and 
CHH estimates. The helium-helium potential energy is, instead, strongly dependent on the wave- function. For the 
uncorrelated cases the repulsive core is overwhelmingly dominant, except in the Pauli suppressed (11)~ state. The 
short-range structure of the correlated and CHH wave functions results in a slightly attractive value of (Vno-He) (from 
-0.3 K to -1.6 K). 

For the mixed ^He-^He-SFg trimer we give only the energies of the lowest lying (0)+ state. The extension of the 
CHH theory presented in Section III to this type of trimer is straightforward. All energies consistently sit in between 
the lighter ■^He2-SFe and the heavier *He2-SFg cases. 

It is worth noticing that the strong supression of the mutual He-He repulsion due to the correlations translates 
into a total binding energy which is very close to the sum of the two dimer energies. The sum of the i — energies 
of Table m gives —54.67, —61.13 and —57.90 K, respectively, for the combinations '^He2, ''He2 and '^He-^'He, which 
are very close to the binding energies of the trimer. The practical effect of the correlations is to reduce the He-He 
interaction to a small attraction of w 0.1 — 0.2 K. This result is very different from the findings of Ref. 14, with total 
binding energies smaller than ours by roughly a factor of two. 

We use a simple argument to estimate the accuracy of the infinite dopant mass approximation. From Table III we 
observe that, to a very good extent, the He2-SF6 trimer can be considered as the superposition of two independent He- 
SFg dimers. Accordingly, the a-/9-trimer corrected kinetic energy results r("''3) = T^^jflil + l/2(m„/Af + mp/M)]. 
This correction corresponds to a modification of the total trimer energy less than 1%. 

Structural properties of the trimers are shown in Table IIVI We give the root-mean-square (rms) dopant-helium 
distance, yJ{E?), the rms helium-helium distance, VV12) j ^'^'^ average value of the cosine between the two D-He 
radii, (cos(6'i2)), in the three approaches. ^/{R?) is not very sensitive to the introduction of the He-He correlation. 
For a given trimer, it assumes essentially the same values in the different states, reflecting the fact that the D-He wave 
functions, (jfis and (jfip, are similar. Since the He-SFg interaction is the same for both helium isotopes, the smaller 
■^He mass would produce a larger kinetic energy than ^He. In order to minimize the energy, this tendency is partially 



compensated by a larger value of y (^d_3Hc) with respect to y (^d_4Ho)- ^^'^ He- He average distance increases in 
going form the uncorrelated to the variational and CHH cases in the spatially symmetric states, as a consequence of 
the intoduction of the correlation, which suppresses short-range He-He configurations. The effect is not very visible in 
the spatially antisymmetric state, (11)~, being these configurations already largely inhibited by the Pauli principle. 
We even find a small decrease of VV12) i'^ ^^i^ state after solving the CHH equations. The values of the uncorrelated 
average He-He cosine, (cos(6'i2)), are immediately understood in terms of the structure of $(^3)7, (1, 2). Correlations 
change these values, mostly in the spatially symmetric 5 = states. For instance, the average cosine corresponds to 
6*12 = 0.39 TT in the uncorrelated (10)" state of the ^He2-SF6 trimer, and to 6*12 = 0.45(0.44) tt for the Jastrow (CHH) 
case. The change for the spatially antisymmetric state, (11)~, is much less evident. 
The one-body helium densities (OBD), 
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normalized as 

''dripW(ri) = 1, (30) 



are shown in Figure O for the L^O and 1 states of the ^He2-SFg and ^He2-SFg trimers, obtained in the uncorrelated 
(symbols) and CHH (lines) approaches. The OBDs are very little affected by both the introduction of the He- 
He correlation and by the optimization of the D-Hc wave function. Actually, they are similar to the dimer radial 
probability densities shown in Fig. ^ As in the dimer case, the "^He atom is more localized than the "^He one because 
of its larger mass. For a given trimer, the OBDs do not appreciably depend on the i-values, since the Is and Ip D-He 
wave functions are almost coincident. 

Differences between the various approaches show up in the helium- helium two-body density (TBD), defined as: 

In Figure 21 we display the center-of-mass integrated TBD, 

(ri2) = / dRi2p(2Hri,r2) (32) 



normalized as 

'"dri2p^')(ri2) = f , (33) 

where R12 = ("iiri -I- m2r2)/(mi -I- 7712) is the He- He center-of-mass coordinate and ri2 = ri — r2 is the He-He 

(2) 

distance. p\ {ru) gives the probability of the two helium atoms being at a distance ri2 apart. 

The Pauli repulsion suppresses p^\ri2) in the uncorrelated (11)^ state of the ^He-SFg trimer at short He-He 
distances, in contrast with the uncorrelated (00)"'' one. This behavior makes the former state bound and the latter 
unbound. We recall that the repulsive core of the He-He interaction is Rc ^ 2.5 A. The introduction of the He-He 
correlation depletes the TBD at small r-values in all the states, which result, as a consequence, all bound. In both 
trimers the helium atoms are more closely packed in the spatially symmetric, L — 1 states, consistently with the 
values of ^/(j^fT)' shown in Table Hvl As expected, the spatially antisymmetric (11)^ state is the most diffuse. The 
uncorrelated long-range structures of the TBDs remain essentialy untouched by the correlations. 

We finally show in Figure |5l the TBD around the SFg molecule in the isosceles configurations, p'^\ri,r2 = 7'i,?'i2), 
for the uncorrelated and CHH (00)+, (10)~ and (11)~ states of ^He2-SF6. All of the TBDs vanish at low ri values 
because of the strong D-He repulsion. The isotropic distribution shown by the uncorrelated (00)'*' TBD disappears 
after introducing the He-He correlation, which suppresses the density at low inter-helium distances. As already noticed 
in Fig. 0] the (11)~ TBD is the least affected by the correlations since it displays a short range He- He repulsion due 
to the Pauli principle. 



V. SUMMARY AND CONCLUSIONS 



The study of He2-SF6 trimers reveals the crucial role of the dopant heavy molecule in binding these systems. In 
fact, the SFg molecule acts as a fixed center of force in which the He atoms are moving. We have shown that the 
specific features of the SFg-He interaction gives a rotational band in the excitation spectrum of the dimers. The 
different masses of the '^He and '*He explain in a qualitative way the particular features of each dimer. The solution 
of the dimer Schrodinger equation provides the single-particle wave functions needed to build the trial wave function 
to be used in the variational study of the trimers. 

The ground state energies of the three isotopic trimers, namely ''He2-SF6, ^He2-SF6 and '^He-^He-SFg, have been 
estimated employing a Jastrow correlated wave-function built up as the product of the dimer wave functions times 
a two-body correlation function of the McMillan type between the He atoms, having a single variational parameter. 
The accuracy of the variational approach has been tested against the Correlated Hyperspherical Harmonics expansion 
method. We have found that the variational results are in excellent agreement with the CHH ones at convergence. 

The role of the Jastrow correlation is crucial in order to overcome the strong repulsion between the He atoms. 
Actually the uncorrelated variational approach docs not bind the trimers, except the ■^He2-SF6 one in the (11)~ 
configuration. The reason is that its wave function is spatially antisymmetric, and therefore the two "^He atoms are 
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kept already apart by the Pauli repulsion. Wc stress that the preferred spatial configuration assumed by the two He 
atoms is such that they take advantage from the mutual attraction, suppressing, as much as possible, the short range 
repulsion. As a result, the optimal configuration is not a linear one, with the SFg molecule in the middle of the two 
He atoms. Instead, the He atoms are closer, and their position vectors with respect to the SFe molecule form an 
angle between ~ 70° and ~ 86°, depending on the particular state. The practical effect of the correlations is that the 
binding energy of the trimer is slighlty larger than the sum of the hinging energies of the corresponding dimers. 

The He-He correlation does not particularly affect the the helium probability densities in the trimers, which are 
similar to those in the dimer. In contrast, the correlation is essential in inverting the energy hierarchy between the 
spatially symmetric and antisymmetric configurations. In fact, in the correlated ■^He2-SF6 trimer the (00)^ state, 
symmetric in space with both He atoms in the Is state and 5 = 0, is more bound than the (11)~ state, antisymmetric 
in space with one atom in Is and the other in the Ip state and S = 1 (aligned spins). The uncorrelated approach 
does not even bind the (00)"'" trimer, whereas the (11)~ one is still bound. 

The good agreement between the variational and the CHH estimates for the trimers makes us confident that medium 
size He doped clusters can be accurately described by means of a correlated variational wave function. This trial wave 
function would be built up from the single-particle wave functions obtained after solving the dimer case, and from an 
appropriate Jastrow factor to properly take into account He-He correlations. In this respect. Variational Monte Carlo 
and Fermi Hypernetted Chain techniques seem to be the most likely candidates to microscopically address the study 
of medium-heavy doped helium nanodroplets. 
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-27.336 


-39.184 


4.622 


4.646 


-1.940 


-30.563 


-41.509 


4.549 


4.568 


-4.292 


1 


-26.561 


-39.104 


4.626 


4.650 


-1.541 


-29.964 


-41.459 


4.552 


4.571 


-3.920 


2 


-25.016 


-38.939 


4.635 


4.659 


-0.775 


-28.767 


-41.355 


4.557 


4.577 


-3.184 


3 


-22.707 


-38.677 


4.648 


4.673 




-26.976 


-41.193 


4.566 


4.585 


-2.105 


4 


-19.647 


-38.298 


4.667 


4.693 




-24.598 


-40.963 


4.577 


4.597 


-0.721 


5 


-15.854 


-37.771 


4.693 


4.720 




-21.641 


-40.653 


4.592 


4.613 




6 


-11.354 


-37.047 


4.728 


4.758 




-18.117 


-40.243 


4.612 


4.634 




7 


-6.183 


-36.041 


4.776 


4.809 




-14.040 


-39.707 


4.637 


4.660 




8 


-0.394 


-34.579 


4.848 


4.887 




-9.431 


-39.001 


4.669 


4.693 




9 












-4.318 


-38.058 


4.711 


4.738 





TABLE I: Observables of the bound states of the He-SFe dimers. Energies are in K and distances in A. 





(00)^ ground state 


(00)^ first excited state 




1 


2 


3 


4 


1 


2 


3 





-54.269 


-54.652 


-54.656 


-54.656 








1 


-54.381 


-54.743 


-54.746 


-54.746 


-51.652 


-51.674 


-51.674 


2 


-54.397 


-54.777 


-54.780 


-54.780 


-52.384 


-52.416 


-52.416 


3 


-54.413 


-54.789 


-54.792 


-54.792 


-52.411 


-52.439 


-52.439 


4 


-54.414 


-54.792 


-54.795 


-54.795 


-52.425 


-52.454 


-52.454 


5 


-54.416 


-54.793 


-54.796 


-54.796 


-52.426 


-52.455 


-52.455 


6 


-54.416 


-54.793 


-54.796 


-54.796 


-52.427 


-52.456 


-52.456 


7 


-54.417 


-54.793 


-54.796 


-54.797 


-52.427 


-52.457 


-52.457 


8 


-54.417 


-54.793 


-54.796 


-54.797 


-52.427 


-52.457 


-52.457 



TABLE II: Convergence pattern of Ef for the two lowest lying states of the (00)^ band in ■^He2-SF6 as a function of the 
numbers {n\,nz) of Jacobi polynomials in the expansion I25II for the amplitude (1,2) and 3, respectively. The number of 
Laguerre polynomials is independently optimized for each case. The variational parameter /3 is fixed to 1.59 A. Energies in K. 
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^He-^Hc 


*He- 


^He 


3He-*He 




(00)+ 


(10)- 


(11)- 


(0)+ 


(1)- 


(0)+ 


b 


1.16 


1.17 


1.06 


1.16 


1.17 


1.16 




837.18 


1704.7 


-33.378 


931.84 


1899.2 


881.18 


J^(var) 


-54.746 


-53.859 


-54.078 


-61.346 


-60.783 


-58.045 


j^(CHH) 


-54.797 


-53.871 


-54.161 


-61.442 


-60.843 


-58.102 




23.696 


24.390 


24.390 


21.892 


22.440 


22.794 


rp(var) 


23.936 


24.948 


24.389 


22.074 


22.869 


23.004 


rp(CHH) 


24.301 


25.587 


24.553 


22.668 


23.713 


23.528 


, ,(unc) 
^ D-Ha 


-78.368 


-78.288 


-78.288 


-83.018 


-82.967 


-80.693 


■,r{var) 
^ D-He 

,r(CHH) 


-78.045 
-78.392 


-77.661 
-78.252 


-78.213 
-78.402 


-82.761 
-83.058 


-82.469 
-82.988 


-80.401 
-80.702 


"He-He 

"He-He 
y(CHH) 
"He-He 


891.85 
-0.638 
-0.705 


1758.6 
-1.146 
-1.206 


20.520 
-0.254 
-0.312 


992.97 
-0.660 
-1.052 


1959.7 
-1.183 
-1.569 


939.08 
-0.649 
-0.908 



TABLE III: Total, kinetic and potential energies (in K) of the He2-SF6 trimers in the uncorrelated, variational and CHH 
approaches. 6 is the adimensional parameter of the McMillan correlation function. 





3 


Hc--''He 


*Hc- 


^Hc 


'^He-^He 




(00)+ 


(10)- 


(11)- 


(0)+ 


(1)- 


(0)+ 


r- — — ( unc) 


4.646 


4.648 


4.648 


4.568 


4.570 


4.646/4.568 


n TrriCHH) 


4.655 

4.645 


4.665 

4.645 


4.650 

4.653 


4.576 

4.568 


4.584 

4.569 


4.656/4.575 

4.646/4.567 




6.570 


5.382 


7.580 


6.461 


5.288 


6.516 


V(^?2> 


6.996 
6.845 


6.029 
5.903 


7.642 
7.476 


6.893 
6.351 


5.946 
5.363 


6.944 
6.533 


(C0S(6»12))<""'=' 





1/3 


-1/3 





1/3 





(C0s(ei2))('^'°'') 


-0.133 


0.163 


-0.355 


-0.138 


0.157 


-0.135 


(C0s(ei2))(^^^' 


-0.087 


0.195 


-0.298 


0.033 


0.260 


-0.006 



TABLE IV: Root mean square dopant-helium, and helium-helium, ^J{rl^, radii and average cosine between the two 

dopant-helium radii, (cos(6i2) in the uncorrelated, variational and CHH approaches. In the first two rows of the last column, 
the first (second) value gives the D-^He (D-^'He) rms radius. Distances in A. 
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FIG. 1: Dimer radial probability densities (^He-SFe: upper panel; ^He-SFe: middle panel); He-SFg interaction (lower panel). 
The probability densities correspond to the quantum numbers £ = (solid line), £ = 4 (dashed line) and £ = 8 and 9 (dash-dotted 
line) for ^He and ''He, respectively. Distances in A, potential in K and densities in A~^. 



14 



-1 I I I I I I I r 



He-SF, / 
" / 
/ 

/ 

/ 

/ »■ 

/ .•■ He-SF, 
/ • '' 

/ 

/ 

/ 

/ .* 



/ 

/ 

/ 

/ «• 

/ . 
/ •■ 



I 

/ : 
/.* 
/.' 

/■* 



:)£ I I I I I I I L_ 

20 40 60 80 

1(1+1) 



100 



FIG. 2: Differences eo; — eoo (in K) vs. £{£ + 1) for the the He-SFe dimers. Squares correspond to the ^He-SFe results and 
stars to the '^He-SFe ones. The Unes are Unear fits to the numerical values. 
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FIG. 3: Helium densities in the L—0 (upper panel) and L—1 (lower panel) states of the He2-SF6 trimcrs. The symbols 
denote the uncorrelated densities, the lines stand for the CHH ones. In the upper panel the circles and the dashed line stand 
for ^He2-SF6 while the triangles and the dash-dotted line correspond to ■^He2-SF6. The same notation is used for the L=l 
*He2-SF6 case in the lower panel. For the L=l ^He2-SF6 trimer two states are reported, corresponding to S=0 (triangles and 
dotted line) and S=l (stars and dash-dotted line). However, at the scale of the Figure they are hardly distinguishible. 
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FIG. 4: Helium-helium center-of-mass integrated two-body densities in the ^He2-SF6 (upper panel) and *He2-SF6 (lower panel) 
trimers in the uncorrelated and CHH approaches. The uncorrelated (10)~ two body density for ^He2-SF6 is very close to the 
(00)"*" one and it is not shown in the Figure. The same holds for the uncorrelated (1)~ TBD for *He2-SF6. 
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FIG. 5: Helium-helium two-body densities in the isosceles configurations, p'^\r\,r2 — ri,ri2), for the ^Hc2-SF6 trimer in the 
uncorrelated (00)+ (a), (10)" (c) and (11)" (e) states, and in the CHH (00)+ (6), (10)" (d) and (11)" (/) states. Distances 
in A and densities in A~^. 



